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1 .  INTRODUCTION 

This  memorandum  describes  a  method  for  computing  the 
solution  to  the  boundary  value  problem  consisting  of  a  plane 
sound  wave  incident  on  an  infinite-length  right  circular 
cylinder  and  a  section  of  a  right  circular  cylindrical  shell 
concentric  with  the  cylinder.  The  problem  is  formulated  so 
that  the  section  of  the  cylindrical  shell  shields  the  circular 
cylinder  from  the  incident  plane  wave.  The  frequency  of  the 
incident  wave,  the  surface  impedance  of  the  cylinder  and  shell, 
the  relative  spacing  of  the  cylinder  and  shell  and  the  angular 
extent  of  the  shell  are  parameters  of  the  solution.  This 
boundary  value  problem  has  been  considered  because  it  repre¬ 
sents  a  tractable  formulation  for  investigating  the  sound  field 
in  the  vicinity  of  a  sonar  transducer,  produced  by  a  source, 
such  as  the  ship's  screws,  after  this  sound  field  has  been 
modified  by  a  baffle  interposed  between  the  source  and  the 
transducer.  The  cylinder  represents  the  transducer  and  the 
section  of  a  cylindrical  shell  represents  the  baffle. 

The  material  in  this  memorandum  is  concerned  primarily 
with  the  development  of  a  formal  solution  to  the  problem  as 
formulated  above.  Justification  for  applying  this  mathematical 
model  to  the  baffle- transducer  complex  and  a  determination  of 
the  range  of  parameters  which  can  be  usefully  considered  will 
be  considered  in  a  subsequent  memorandum. 
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II .  FORMULATION  OF  THE  PROBLEM 

Consider  the  problem,  shown  schematically  in  Figure  1,  of 
a  plane  sound  wave  incident  on  a  right  circular  cylinder  of  in¬ 
finite  '  length  and  radius  a,  and  a  section  of  a  right  circular 
cylindrical  shell,  also  of  infinite  length,  located  at  radius  b, 
where  b  >  a.  The  shell,  which  will  be  called  the  baffle  in  the 
following,  has  an  angular  width  2  (tt  -  cp)  and  is  positioned  with 
respect  to  the  incident  wave  so  that  symmetry  in  the  angle  8  is 
preserved.  (The  case  in  which  symmetry  in  the  angle  8  is  not 
preserved  is  discussed  in  Section  V.)  The  thickness  of  the 
baffle  is  assumed  to  be  negligible.  It  is  desired  to  obtain  a 
solution  to  the  wave  equation  which  describes  the  total  sound 
field  throughout  all  of  space  external  to  the  cylinder  for  the 
incident  plane  wave  and  the  scattered  and  diffracted  waves 
arising  from  specified  boundary  conditions  on  the  surface  of 
the  cylinder  at  radius  a  and  on  the  inner  and  outer  surface  of 
the  baffle  at  radius  b. 

As  is  well  known,  the  wave  equation  in  cylindrical  co¬ 
ordinates  for  cases  in  which  the  solution  is  not  a  function  of 
axial  position*  is 


+  __1  =  1 

dr  2  36  2 

r  c 


(1) 


where  p  is  the  sound  pressure,  r  and  0  are  the  ihdependent  co¬ 
ordinates,  t  is  the  time  and  c  is  the  speed  of  sound.  The 
general  solution  to  this  equation,  which  can  be  obtained  by 
assuming  that  the  solution  has  the  form  of  a  product  of  functions 
of  each  of  the  independent  coordinates  and  the  time,  is,  for  a 

*In  this  boundary  value~probiem  the  baffle  and  cylinder  are  as¬ 
sumed  to  be  of  infinite  length  in  order  that  the  end  effects 
may  be  neglected.  The  solution  is  therefore  independent  of  axial 
position. 
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Fig.  I -GEOMETRY  OF  THE  INFINITE-LENGTH 
CYLINDER  AND  CYLINDRICAL  SHELL 
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sinusoidal  time  dependence, 


p  =  V'  [A  J  (kr)  +  B  N  (kr)]  cos  n9  e  (2) 

/  ,  nn'  nn 

n=0 

where  An  and  are  undetermined  constants,  Jn(kr)  and  Nn(kr) 
are  the  Bessel  and  Neumann  functions  respectively,  (i.e.,  the 
two  independent  solutions  to  Bessel's  equation  for  integral 
orders  ),  u)  is  the  angular  frequency  of  the  incident  wave, 
k  =  2tt/X  is  the  wave  number  and  \  is  the  wavelength  of  the  in¬ 
cident  wave.  In  order  to  obtain  a  better  physical  insight  in¬ 
to  the  meaning  of  this  solution,  it  can  be  rewritten  in  the 
f  orm^ 

CO 

P  =  £  [CnHn(1)(kr)  +  DnHn(2)  (kr)l  cos  n9  e_:LU)t,  (3) 

n=0 


where  H  ^^(kr)  =  J  (kr)  +  i  N  (kr)  is  the  Hankel  function  of 
the  first  kind  and  represents  a  cylindrical  wave  propagating  radi¬ 
ally  outward;  and  Hn^2^(kr)  =  J^(kr)  -  i  Nn(kr)  is  the  Hankel  func¬ 
tion  of  the  second  kind  and  represents  a  cylindrical  wave 

3 

propagating  radially  inward  ,  corresponding  to  reflected  and 
incident  cylindrical  waves  respectively.  The  radial  particle 
velocity  is  obtained  from  the  pressure  according  to  the  relation 


^See  e.g.,  P.  M.  Morse,  Vibration  and  Sound,  McGraw-Hill,  New 
York,  1948,  p.  298. 

2 Ibid,  pp.  188,  196. 

3 

P.  M.  Morse  and  H.  Feshbach,  Methods  of  Theoretical  Physics, 
McGraw-Hill,  New  York,  1953,  p.  1371. 


3 


TRAC  OR,  INC. 


1701  Guadalupe  St. 


Austin  I,  Texas 


ipck  v  =  ,  (4) 

where  v  is  the  radial  particle  velocity,  i  =^/-l,  p  is  the  density 
of  the  medium  and  the  remaining  terms  have  been  previously  de- 
f ine  d . 

Before  proceeding  to  an  evaluation  of  the  undetermined 
constants  so  that  the  specified  boundary  conditions  are  satisfied, 
it  is  of  interest  to  recall  briefly  several  fundamental  identities 
and  to  discuss  some  notation  pertinent  to  representing  boundary 
conditions  on  the  baffle.  Consider  a  function  f (0 )  in  the  form 
of  equations  (2)  and  (3) . 


f (0)  An(r)  cos  n0> 

n=0 


o<|  e  |  <cp 


(5) 


=  0 


cp<  |  0  |  <rr 


The  function  f(8)  can  be  written  as  a  single  expression  valid 
for  all  0  by  expanding  it  in  a  Fourier  series.  The  Fourier 
series  expansion  of  f(0)  results  in  the  identity 

CX>  CO 

f(9)  =  2ttcp  ^  Aj  (r){^j-c°s  n0 ,  0<  1 0 1  <tt  (6) 

n=0  j=0 

where  €q  =  1,  6  =  2  when  n  ^  0,  and 

/  n\=  sin  |n  -  11  cp  ,  sin  (n  ±  i )  cp 
UJ  I  n  -  j  cp  (  n  +  j)  cp 
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Here  the  Fourier  coefficients  have  been  determined  in  the  usual 
4  . 

manner,  i.e.. 


^  I  f (0 )  cos  n9  d8 


which  in  this  case  reduces  to 


Cp  ra 


^  0  X  Aj^  C°S  C°S  n6  de  =  2n  ^  X  Aj^  {j} 


As  indicated  in  eq .  (6)  this  expansion  is  valid  for  all  0. 

The  expression  for  f (6 )  given  in  eq.  (6)  may  also  be  ob¬ 
tained  as  the  product 

f(0)  =  g<9)h(6),  (7) 


where 


g(9)  =  1, 


o<  |  e  |  <cp 


cp<  I  6  |  <TT 


h(6)  -  £  An<r)  cos  n0, 


o<  1  e  I  <tt  . 


See  e.g,,  D.  Jackson,  Fourier  Series  and  Orthogonal,  .Polynomials 
Mathematical  Association  of  America,  1941,  pT  7. 
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h<e)  -  X  Bn(r)  cos  n9  ,  0<|9|<rr. 

n=0 

It  is  convenient  to  develop  the  solution  for  two  different 
geometric  regions.  Let  Region  1  be  defined  by  a<r<b,  0<|b|<tt, 
and  let  Region  2  be  defined  by  r>b,  0< |  9  |  <tt  .  In  Region  2  the 
pressure  p2  can  be  written  as  the  sum  of  the  incident' plane  wave 
plus  a  general  set  of  cylindrical  waves  traveling  radially  out¬ 
ward.  Thus 


P2  (r ,  9 )  =  £€ninjn(kr)  +  (kr)  ^cos  n0  e 

n=0  (12) 


where  the  first  term  under  the  summation  is  the  expression,  in 
cylindrical  coordinates,  of  a  plane  wave^  traveling  as  shown  in 
Figure  1,  €  =  1,  =  2  when  n  ^  0,  and  i  =  V^T.  In  Region  1 

the  pressure  p-^  can  be  written  as  the  sum  of  a  general  set  of 
cylindrical  waves  which  are  propagating  radially  inward  and  a  set 
which  are  propagating  radially  outward  as  given  by  eq .  (3). 

Thus 


CD 

pl<r>0)  =  Z  [BnHn(1)(kr)  +  CnHn(2)(kr)  cos  n0 


-iuut 


n=0 


(13) 


At  the  boundary  r  =  a,  i.e.  on  the  surface  of  the  cylinder,  either 
p-^(a,0)  or  the  radial  particle  (or  medium)  velocity  v^(a,0)  is 
a  prescribed  function.  On  the  closed  portion  of  the  boundary  on 
the  baffle  at  r  =  b,  i.e.,  for  cp<  |  9  j  <tt  ,  either  p^(b,9)  or  v^(b,9) 
is  a  prescribed  functioia,  while  on  the  open  portion  of  the 

“*P.  M.  Morse,  loc.  cit.,  p.  347. 
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boundary  at  r  =  b,  i.e.  for  0<  1 9  |-  <cp ,  both  the  pressure  p1(b,0), 
and  the  radial  particle  velocity  v^(b,0),  must  equal  the  pressure 
P2(b,0),  and  velocity  v2(b,8),  respectively,  in  order  to  es¬ 
tablish  continuity  in  the  pressure  and  radial  particle  velocity 
between  Regions  1  and  2. 

For  the  pressure  and  velocity  in  Region  1  the  boundary 
conditions  at  the  radius  r  =  b  can  be  written  as 

P1  (b ’  0  )  ={pi_(b>e)  ■  ^ Pi  (b  » e  )  D o<  |  g  |  <cp  }  +  tp2  (b,  0)1q<  I  0  |<cp. 

(14) 

Vl(b,0)  ={v1(b,0)  -  [v1(b,8)]()<|  0  j  }  +  Cv2(b,0)]o<|0|<cp 

(15) 

where  the  brackets,  [],  denote  a  Fourier  expansion  as  given  by 
eq .  (6),  and  the  quantities  without  brackets  are  applicable  for 
all  0.  Braces, {}  ,  are  used  to  indicate  the  terms  which  repre¬ 
sent  a  Fourier  expansion  as  given  by  eq .  (10).  The  corresponding 
equations  at  the  radius  r  =  b  for  the  pressure  and  velocity  in 
Region  2  are 


P2(b,0)  =  {  p2(b,0) 

[P2(b>0)^o<|  0  |  <cp  . 

S’  +  ^Pi(b»e)^o<j  0  |  <cp 

(16) 

v2(b,0)  =  {  v2(b,0) 

-  t>2(b,e))o<|0|<cp  } 

+ 

1 — i 

< 
h- 1 

"c? 

CD 

1 — I 

o 

A 

CD 

A 

-6 

(17) 

The  constants  A  ,  B  and  C  will  be  evaluated  for  a  specific 
n’  n  n  r 

set  of  boundary  conditions  in  Section  III,  corresponding  to  a 
physical  situation  in  which  the  outside  of  the  baffle  is  perfectly 
rigid,  i.e.,  v2(b,0)  =  0  for  cp<  |  0 1  <tt  ;  the  inside  of  the  baffle 
is  soft,  i.e.,  p^(b,0)  =  0  for  cp<|0|<n;  and  the  cylinder  is  rigid, 
i.e.,  v-^(a,8)  s  0.  In  Section  IV  the  solution  will  be  evaluated 
for  more  general,  non-zero  values  at  the  boundaries. 
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III.  EVALUATION  OF  THE  UNDETERMINED  CONSTANTS 

The  undetermined  constants,  A  ,  B  ,  and  C  must  be  evalu- 

’  n’  n  n 

ated  in  order  to  satisfy  the  boundary  conditions  specified  on 
the  cylinder  and  both  faces  of  the  baffle.  The  case  to  be  con¬ 
sidered  in  this  section  is  that  of  a  rigid  outer  baffle  surface, 
a  soft  inner  baffle  surface,  and  a  rigid  cylinder  surface.  A 
formal  statement  of  these  boundary  conditions  is: 


v2 (r , 0 ) 


=  v1(r,0) 
r=b  1 


y 

r=b 


0<  |  9  |  <cp  (18) 


=  0  , 


Cp<  1  0  l  <TT 


Px  (r ,  0 ) 


r=b 


P2(r»e) 


r=b 


0<  |  8  |  <cp  (19) 


=  0  , 


cp<  |  6  |  <rr 


(r , 0  ) 


r=a 


0  , 


0<  |  9  |  <rr  (20) 


The  solutions  for  Regions  1  and  2  can  be  separated  into  Bessel 
and  Neumann  functions  .  Equations  (12)  and  (13')*  then  become 


If 


£  inJ  (kr)  +  A  J  (kr) 
n  n  ■  '  n  n '  / 


+  iAnNn(kr) 


cos  n0 


(21) 


n 


=0 


?The  time  component  e  *'u,t  of  the  solution  will  henceforth  be 
suppressed. 
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and 

CO 

P-L  =  [(Cn  +  Dn)  dn(kr)  +  i(Cp  -  Dr)  Nn(kr)]  cos  n6  (22) 

n=0 

Without  loss  of  generality,  eq .  (22)  can  be  rewritten  as 

CO 

Pj_  =  ^  C En Jn ^kr ^  +  iFnNn(kr)^  cos  n0  (23) 

n=0 

where  C  +  D  =  E  and  (C  -  D  )  =  F  . 
n  n  n  v  n  n'  n 


Consider  the  boundary  condition  given  in  eq .  (20).  The 
radial  particle  velocity  for  Region  1  is  obtain  by  using  eq .  (4). 


ipckv-^ 


r=a 


k[ (ka )  +  iFnN^(ka)]  cos  (n0) 

n-0 


0 


(24) 


where  the  primes  indicate  differentiation  with  respect  to  the 
argument  kr  and  the  Bessel  and  Neumann  functions  have  been  evalu¬ 
ated  at  r  =  a.  Equation  (24)  implies  that 


J' (ka) 

F  =  i  E  — -  . 

n  n  N;(ka) 


Thus,  the  solution  in  Region  1  can  be  written  as 


(25) 


Pi 


n=0 


Jn(kr) 


J^Cka) 

(ka) 


Nn(kr) 


cos  n8  . 


(26) 


Consider  now  the  substitution  of  the  boundary  conditions, 
equations  (18)  and  (19),  into  the  identities  given  in  equations 
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(14)  and  (17).  These  equations  reduce  to 
P]_(bj6)  =  C  P2  Ct* ,0)3 Q< |  e  j <cp 


(27) 


and 


v2(b,0)  =  [v1(b,e)]o<|0|<cp  (28) 

The  identities  expressed  in  equations  (14)  and  (17)  were  chosen 
in  this  case  because  the  quantities 

{Pl(b,  6)  -  Cpx  (b ,  0  )  3^  |  ©  |  <cp  }  and|v2(b,0)  -  Cv2(b,0)]o<|0|<q3| 

were  specified  as  boundary  conditions.  Equations  (14)  and  (17) 

therefore  lead  to  a  direct  determination  of  A  and  E  .  If,  for 

ti  n 

example,  the  velocity  were  specified  on  the  inside  of  the  baffle 

and  the  pressure  were  specified  on  the  outside  of  the  baffle, 

then  equations  (15)  and  (16)  would  be  more  advantageous  in 

solving  for  the  undetermined  constants  A  and  E  . 

D  n  n 

The  quantities 

[p2(b,e)]0<|e|<cp  and  [v1(b,e)3o<|0|<p 

are  determined  according  to  eq.  (6)  using  the  forms  of  the 
solution  for  Regions  1  and  2  given  in  equations  (21)  and  (26) . 
Equations  (27)  and  (28)  then  become 


n=0 


J 1 (ka) 

Jn(kb)  -  -n-~  -■  N„(kb) 


^n (ka ) 


n 


cos  ne  =  Z  2S 9  Zn  iJjj(kb) 

n=0  j=0 


+  Aj-Jj  (kb)  +  iAj Nj  (kb)j-  V  j-  cos  n0 


(29) 


11 
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A  =  - - - 

n  J^(kb)  +  iN^ (kb) 


-  €  inj' (kb) 
n  n v  ' 


This._value  of  An  is  substituted  into  eq.  (31)  to  yield 


3  J 
n  n 


J' (ka) 

(kb)  -  -a -  N  (kb ) 

N'(ka)  n 


.  >1  a  (kb)  +  iN.  (kb)  €.  ^ 

3J.  (kb)  +  — ■ - ] -  ^r1  cp  ) 

3  J' (kb)  +  IN! (kb)  L  2  ^ 

J  J  t«0 


{Jhkb)  '  KS  Ni(kb)}  0}  • 


The  order  of  summation  of  eq .  (34)  can  be  interchanged  and  the 
equation  rearranged  as  follows: 


J' (ka) 


En  Jn(kb)  "  ~T -  N  (kb)  = 

n  L  n  N^ka)  n  J 


6  y-'  r  Jl(ka)  -j  y— v  /-•->  J .  (kb)  +  IN.  (kb)  <-  ■ 

cp  >  E  J  J I  (kb)  -  — -  N !  (kb)  i  V  JU  -1 - 1 - cp  J  J 

2tt  /Lv  [ _  l  N !  (ka )  1  J  ^  W  j!  (kb)  +  iN!(kb)  2tt  \3. 


N|(ka) 


N;(kb)}  Z{i}^ 

j=0  3 


J,(kb)  +  IN, (kb)  €. 


• (kb)  +  iN! (kb) 


^,CP  j 


€  f  .  .  J.(kb)  +  iN.(kb)-)  r  q 

+  y^cp  V  -U.iJj.(kb)  -  €.iJj!(kb)  -Ll— - 12 -  {?} 

2tt  J  3  3  3  j!(kb)  +  iN.'(kb)  J  UJ 

i=n  J  J 
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The  summations  over  the  index  j  involve  known  constants  only  (the 
Bessel  and  Neumann  functions)  so  that  the  index  j  can  be  summed 
out.  Equation  (35)  therefore  can  be  written  in  the  form 

En  -  (35a) 

l 

which  is  an  infinite  set  of  inhomogeneous  equations  relating  the 
E's  to  each  other.. 

A  qualitative  examination  of  the  functional  forms  of  the 
quantities  and  (3^  shows  that  the  factors  and  {j}  will 

have  a  tendency  to  cause  the  major  contribution  to  the  sum  on  t 
to  arise  for  values  of  t  in  the  vicinity  of  n.  One  would  expect, 
therefore,  that  it  should  be  possible  to  obtain  accurate  numerical 
values  for  the  first  N  of  the  E's  if  the  indices  in  the  summation 
(n  and  -f)  in  eq ,  (35a)  are  restricted  to  a  finite  number  M  >  N. 

If  the  indices  are  restricted  in  this  manner,  eq .  (35a)  becomes 
a  set  of  M  inhomogeneous  equations  in  M  unknowns ,  which  can  be 
evaluated  numerically.  The  pressure  can  then  be  calculated  from 
eq  „  (26).  The  number,  N,  of  accurate  E's  required  is  determined 
by  the  convergence  rate  of  eq.  (26).  The  number  M  of  E's  which 
must  be  included  in  eq.  (35)  to  obtain  accurate  values  for  the 
first  N  of  the  E's  is  determined  by  the  convergence  effects 
of  the  factors  { j  ]•  and  in  computing  the  sums  over  the 

index  l. 

The  qualitative  aspects  of  obtaining  numerical  values  for 
the  pressure  described  in  the  preceding  paragraph  have  been  ex¬ 
amined  quantitatively  and  have  been  found  to  be  substantially 
correct.  The  complete  computational  procedure  as  well  as  some 
numerical  results  will  be  described  in  a  subsequent  memorandum. 
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IV .  GENERAL  BOUNDARY  CONDITIONS 

In  the  previous  section  a  solution  to  the  baffle-cylinder 
sound  field  was  developed  for  the  case  of  zero  radial  particle 
velocity  on  the  outside  of  the  baffle  and  zero  pressure  on  the 
inside  of  the  baffle.  Solutions  for  more  general  boundary  con¬ 
ditions  on  the  baffle  can  also  be  developed  using  a  procedure 
identical  to  that  of  Section  III. 

Expressions  for  the  quantities  p^,  v^ ,  p ^ >  and  V£  on  the 
boundary  r  =  b  are  given  in  equations  (14),  (15),  (16)  and  (17). 

In  these  equations  the  quantities  in  the  braces  ||  represent  the 
specified  boundary  conditions.  Any  combination  of  pressure  and 
radial  particle  velocity  can  be  specified  as  boundary  conditions 
on  the  baffle.  For  example  if  the  pressure  on  the  outside  of  the 
baffle  is  specified  as  some  function  p2(0)  and  the  radial  particle 
velocity  on  the  inside  of  the  baffle  is  specified  as  some  function 
V^(0)  then  equations  (15)  and  (16)  become 


°°  TT 


vL  (b ,  e )  =  ^  ^  v]_(®)  cos  n0  de  cos  n0  +  t  v2  (b ,  9  )  3  o<  |  q 


and 

Po  (b, 0) 


n=0 


CO  TT 


|cp 

(36) 


=  ^  J*  P2(8)  cos  nB  de  cos  n9  +  [P1(b,6)]o<|0  |<cp 


n=0 


(37) 


where  the  series  represent  the  Fourier  expansions  of  V-^(9)  and 
P2(0)  respectively.  Here  equations  (36)  and  (37)  are  analogous 
to  equations  (27)  and  (28)  of  Section  III.  The  procedure  for 
determining  the  unknown  coefficients  An  and  En  is  then  identical 
to  that  of  the  previous  section. 
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The  solution  for  the  problem  using  any  other  combination 
of  boundary  conditions,  e.  g.  pressure  specified  on  both  the 
inside  and  outside  of  the  baffle,  is  obtained  by  choosing  the 
appropriate  pair  of  identities  (equations  (14)  -  (18))  and  pro¬ 
ceeding  as  in  Section  III. 
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V.  PLANE  WAVES  OBLIQUELY  INCIDENT  TO  THE  BAFFLE -CYLINDER 

In  Section  III  a  solution  was  obtained  for  the  boundary 
value  problem  consisting  of  a  plane  wave  normally  incident  on 
the  baffle-cylinder.  (See  Figure  1.)  In  this  section  the 
solution  will  be  obtained  for  the  more  general  case  of  a  plane 
wave  incident  on  the  baffle-cylinder  at  some  angle  as  shown 
in  Figure  2.  The  development  of  the  solution  for  the  oblique 
incidence  case  closely  follows  the  normal  incidence  case 
presented  in  Section  III.  The  expansion  for  the  normally 
incident  wave  in  polar  coordinates  is 

00 

p  =  ^  €ninJn(kr)cos  n9 

n=0 


where  the  symbols  are  defined  in  previous  sections .  By  a 
simple  coordinate  transformation,  the  off-axis  wave  shown  in 
Figure  2  can  be  written  as 

CO 

p  =  >  €  inJ  (kr)  [cos  ntcos  n9  +  sin  n'l'sin  n8] 

L-.  n  n 

n=0 


where  t  is  the  angle  between  the  wave  front  normal  and  the 
0=0  axis.  For  the  off-axis  case,  the  solution  is  no  longer 
symmetric  in  0;  consequently,  both  sine  and  cosine  terms  must 
be  present  in  the  Fourier  expansion  of  the  solution.  Thus 
the  forms  of  the  solution  in  Regions  1  and  2  must  be  written 
as 

P.M.  Morse,  loc.  cit .  p.  3^7 

^See  e.g.  R.  V.  Churchill,  Fourier  Series  and  Boundary  Value 
Problems ,  McGraw-Hill  Company,  New  York,  1941,  Chapter  IV. 
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and 


Pl  "  E  [EnJn(kr)  +  iFnNn(kr)]cos  n0 
n-0 

CO 

+  [KnJn(kr)  +  iLnNn(kr) 3 sin  n0 

n-0 


(38) 


X”1  [€  inJ  (kr)cos  n*!1  +  A  J  (kr)  +  iA  N  (kr)]cos  n9 

n  nv  '  n  nv  '  n  nv 

n-0 

(39) 

CO 

+  V'  [€  inJ  (kr)sin  nt  +  G  J  (kr)  +  iG  N  (kr))sin  nQ . 
l_j  n  nx  n  ns  n  n 

n-0 


The  evaluation  of  the  constants  A  ,  E  ,  F  ,  G  ,  K  ,  and  L 

n’  n  n’  n  n  n 

follow  the  procedure  outlined  in  Section  III. 

Consider  the  boundary  conditions  given  by  equations  (18), 
(19) ,  and  (20) .  Substitution  of  equations  (38)  and  (39)  into 
(18) ,  (19) ,  and  (20)  and  performing  the  operations  indicated 
in  Section  III  yields 


E  M  -  ~  cf 
n  n  2rr 


'  I 

l-  0 


€^i'{/J^(kb)cos  l-ty  +  Aj, J^(kb)  +  iA  N  (kb)j  { -t  j- 

(40) 


K  M 
n  n 


€ 

n 

2tt 


A-0 


1 


C^i'/v'J^(kb)sin  •U1  f  G^J^(kb)  4-  iG,N,(kb)[-  1  n 


(41) 
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€  inJ '  (kb)cos  n\|i  +  A  J '  (kb)  +  iA  N '  (kb) 
n  nx  y  n  nv  '  a  nv 


and 


€ninjA(kb)sin  n+  +  GnJ;(kb)  +  iGnNn<kb) 


Z  {vi}{?}  * 

(42) 


1=0 


-  *  Z  M  ft}  •* 

(43) 


1=0 


where 


ft} 


s  in  1 

-t-n 

cp  sin  1 

+  n; 

»cp 

|  -6-n 

cp  1 

+  n) 

icp 

Jn_l  (ka )  “  Jn+l(ka) 

Mn  ”  Jn<kb>  -  N".ftka)  -  Vkb>  - 

J  (ka)  -  J  ,  i  (ka) 

m;  ■  j;<kb>  -  vftka>'  -  vftk*>  Nb<kb)  ■ 

and  all  other  symbols  are  defined  as  in  previous  sections. 
Equations  (40)  and  (42),  (41)  and  (43)  can  be  combined  to  yield 
equations  Similar  to  eq.  (35),  i.e. 


€  00 

EnMn  *  2?  'P  Z  {^cos  ^CJ^(kb)  -  R^J|(kb)]} 


1=0 


(44) 


03  CO 


+  fe'pZ  Z  ’{"}{}}  RjMib  ■ 

j  =0  1=0 
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co 

Knf,n  -  fe  ”  I  {5}  {v**ln  t«Jt(kb)  -  V;<kb)]} 
1=0 


(45) 


CO  05 

€  t  .  . — ,  r  'i  *  r i  *  €  . 


+  5W  v  X  X  {"}  {j}  5*  *  RjMjKj  • 


j  =0  1=0 


where 


R. 

3 


Ji  (kb)  +  iN. (kb) 
(kb)  +  iNj  (kb) 


Equations  (44)  and  (45)  can  be  solved  in  the  same  manner  as 
eq.  (35)  to  yield  values  of  and  Kn  to  complete  the  solution 

for  p^(r,0).  For  boundary  conditions  other  than  those  given 
in  equations  (18) ,  (19) ,  and  (20)  the  procedure  for  the 
treatment  of  more  general  boundary  conditions  given  in  Section 
IV  is  employed. 

The  solution  presented  in  this  section  will  reduce  to 
the  normal  incidence  case  presented  in  Section  III  if  the 
angle  I1  =  0  is  substituted  into  equations  (44)  and  (45)  . 
Equation  (44)  then  reduces  to  eq .  (35)  of  Section  III.  Equa¬ 
tion  (45)  reduces  to  an  infinite  set  of  homogeneous  equations 
in  the  unknown  which  is  satisfied  only  if  Kn=0,  resulting 
in  a  solution  identical  to  that  of  Section  III. 

Complex  incident  sound  fields  can  be  considered  using 
the  solution  presented  in  this  section.  Assume  several  plane 
waves  are  incident  on  the  baffle-cylinder  at  angles  of  in¬ 
cidence  The  scalar  wave  equation  (eq .  1)  is  a  linear  equa¬ 

tion;  hence,  any  linear  combination  of  solutions  to  the  wave 
equation  is  itself  a  solution.  Therefore  the  solution  to  the 
complex  sound  field  containing  M  different  incident  plane  waves 
can  be  written  as 
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Pi  = 


M  “  r 

I  I  {[ 

m=l  n=0 


E  J  (kr)  +  iF  N  (kr) 
mn  n  '  mn  n 


cos  n9 


(46) 


K  J  (kr)  +  iL  NXT(kr) 
mn  nN  mn  N'  ' 


sin  n0 


M  « 


^  Y  a  €  inJ  (kr)cos  nty  +  A  J  (kr)  +  iA  N  (kr)  cos  n0 
/  ,  / ,  |  m  n  n  m  mn  n  mn  n  j 


m= 


=1  n=0  L 


(*7) 


>€ninjn(kr)sin  ntm  +  GmnJn(kr)  +  iGmnNn(kr) 


sin  n0 


where  am  is  an  amplitude  weighting  function  for  the  various 
plane  waves.  Here  it  has  been  assumed  that  the  various  in¬ 
cident  waves  have  the  same  angular  frequency.  However  the 
solution  given  in  equations  (46)  and  (47)  could  easily  be  gen¬ 
eralized  to  include  incident  plane  waves  of  different  frequencies. 
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